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Abstract 

In part 1 (Gouin, jl3j), we proposed a model of dynamics of wet- 
ting for slow movements near a contact line formed at the interface 
of two immiscible fluids and a solid when viscous dissipation remains 
bounded. The contact line is not a material line and a Young-Dupre 
equation for the apparent dynamic contact angle taking into account 
the line celerity was proposed. 

In this paper we consider a form of the interfacial energy of a solid 
surface in which many small oscillations are superposed on a slowly 
varying function. For a capillary tube, a scaling analysis of the micro- 
scopic law associated with the Young-Dupre dynamic equation yields 
a macroscopic equation for the motion of the contact line. The value 
of the deduced apparent dynamic contact angle yields for the average 
response of the line motion a phenomenon akin to the stick-slip mo- 
tion of the contact line on the solid wall. The contact angle hysteresis 
phenomenon and the modelling of experimentally well-known results 
expressing the dependence of the apparent dynamic contact angle on 
the celerity of the line are obtained. Furthermore, a qualitative ex- 
planation of the maximum speed of wetting (and dewetting) can be 
given. 

Key words: contact angle, contact line, hysteresis. 

1 Introduction 

In [13], henceforth referred to as part 1, we proposed a model of non- 
Newtonian fluids for slow movements in the immediate vicinity of the contact 
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line that is formed at the interface of two immiscible fluids and a solid. Fluid 
interfaces are modelled by differentiable manifolds endowed with constant 
capillary energy. Solid surfaces are also regarded as differentiable manifolds 
endowed with a position dependant surface energy (we only consider the case 
of surfaces without surfactant). The kinematics of slow isothermal move- 
ments close to the contact line, revisited in the framework of continuum 
mechanics, required the adherence condition to be relaxed at the contact 
line. The velocity field is discontinuous at the contact line and generates a 
concept of line friction but viscous dissipation remains bounded. 
Simple observations associated with the motion of two fluids in contact with 
a solid wall reveal the following behaviour: depending upon whether the 
fluid on the wall advances or retreats, a variable contact angle is observed. 
The value of this apparent dynamic contact angle, also called simply Young 's 
angle, depends upon the contact line celerity. These observations resolve 
distances which are not shorter than a few microns. 

The most notable unanswered problem is connected with the equation 
controlling the macroscopic motion of the contact line and the justification 
of contact-angle hysteresis when the fluid is advancing or receding on a solid 
wall. 

To verify the accuracy of the model proposed in part 1 and, particularly to 
justify the Young-Dupre equation for the apparent dynamic contact angle 
associated with the notion of line friction, a simple academic example is 
considered for which the equations of motion together with the boundary 
conditions can be integrated in a suitable approximation. This is the case 
of a thin cylinder containing an incompressible fluid separated from air by a 
meniscus. The apparatus is a capillary test-tube with rotation of symmetry. 
In our example, we assume that the solid wall of the tube is endowed with a 
fluid-solid surface energy which oscillates periodically with small variations 
with very short wave length relative to the length of the tube0. The fact 
that we have chosen oscillations superposed on a slowly varying function may 
be justified by the possibility to solve the equations of motion analytically. 
More general forms could also be considered. 

For suitable dimensionless numbers corresponding to slow movements, the 
equations of motion of the liquid and conditions on interfaces take simplified 
forms. The two-fluid interface can then be modelled as a spherical cap. It 

'^Surface energies of solid walls with many small wiggles, arising from small-scale mi- 
crostructural changes, appear often in scientific problems; for example, phase transforma- 
tions, protein folding and friction problems (Abeyaratne, Chu and James, 2 ). The scale 
of microstructural changes is a microscopic one and consequently is of an order smaller 
than the length of the tube. 
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turns out that the microscopic motion of the contact system is governed by 
a simple differential equation. This equation is analytically decoupled from 
those of the liquid motion. When the irregularities of the liquid-wall-surface 
energy vary over a length that is vanishingly small relative to the size of the 
capillary tube, the solution of the microscopic motion tends to a limit which 
is the solution of a new differential equation. This differential equation of 
the macroscopic or homogenized motion is not the limit of the microscopic 
equation. 

The liquid in the capillary tube is controlled by a piston. We deduce the 
motion of the two-fluid interface and we study the behaviour of the apparent 
dynamic contact angle for the advance and retreat of the contact line. In so 
doing a hysteresis phenomenon appears. The deduced results are compared 
with those obtained in the literature with methods of statistical physics and 
experimental measurements. 

2 The capillary tube apparatus 

The apparatus is a cylindrical tube of radius a with a vertical symmetry 
axis Ok. A liquid of volume Vq fills the cylinder above a position deter- 
mined by a piston. In accordance with the hypotheses and notations of part 
1, the liquid is in contact with air through an interface T,2u the constant 
air-liquid-surface energy is (T2 = cab- In the motion, air is considered as 
incompressible. The wall of the cylinder (piston included) is denoted by Ei^. 
The wall is inhomogeneous and the surface energy ui = o"^5 — o'bs (differ- 
ence between the superficial energies of solid-liquid and solid-air) depends 
on the geometrical position on Ei^. The value ai is assumed to be rotation- 
ally invariant. The contact line is the curve Tf connecting the two interfaces 
T,it and T,2t- Due to the axi-symmetric geometry, its representation in fig. 
1 is a point P of which the position is given by the abscissa z; the point J 
of the piston, the position of which is given by the abscissa L, is a function 
of time commanded by an operator. 

On the interval JI of the wall, the surface energy ai is assumed to be 
constant with the value aio- We consider the case when the surface energy 
(Ji on the interval IP of the wall is such that 



where < e ^ 1 is a small dimensionless parameter, | A; |< 1 and kaio is 
positive. 

Let us note that for sufficiently small e, the average value of ui on any 



= <Tio {1 — k sin — ) 
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Figure 1: Circular cylinder with vertical axis filled with a liquid of volume 
Vq held in position by a piston having cap form. The part of the cylinder 
wall with the constant surface energy uio is colored in grey; the part of the 
cylinder wall with the non-uniform surface energy ai is hachured (see main 
text). Referred to the orthonormal axis Oik, z denotes the position of the 
contact line, L = L{t) determines the position of the piston (z = + L). 
The meniscus at the upper free surface is rotational symmetric; the apparent 
dynamic contact angle between the meniscus and the wall is denoted by 6; 
h is the height of the meniscus. 

unspecified macroscopic set of the wall tends to cio when e tends to zero. 
In relation (1) the interfacial energy cJi is built by two terms. The second 

is /e = —k (Tio sin (^z/{ea)^ . The distribution converges to zero when e 

tends to zero and its virtual work on S^i is nul@. We understand that the 

^ Let us calculate the value of the distribution — sin {z/e) when e 0, (e > 0). For 
any function (j> £ C°°[7?, 7?] with compact support, 

< ge,(j) >= sm{x/e)(f){x)dx = e s\nv (j){£v) dv — —e [cos u <;/>(£?;) + 

cosv ((>' {ev) dv. 

The first term [cos v (;/)(ev)]l^ = and |< >| < e'^ f^^ \ <^'(e«) | dv = e / +^ 

(fi'ix) I dx . From fl^°° | (jj'i^) I dx < oo, we obtain hm < g^jfj) > — and lim g^ — 

°° e^0+ £-.0+ 

in a weak sense. 
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associated force has an action only on the contact hne Fj. 
It is proposed to study the motion of the contact hnc, i.e. to determine the 
value of the position z{t) as a function of time. To this end, it is necessary 
to study the liquid motion. 

3 Equations of motion and boundary conditions 

In part 1, we proposed the equations of motion and boundary conditions of 
a capillary motion of two fluids in contact with a solid wall. The equations 
of motion are 



where (/) denotes the volumetric force, p the density, a the acceleration vec- 
tor, Q the (non-Newtonian) viscous stress tensor and p the pressure. At the 
solid wall, adherence conditions are required except at the contact line. On 
the meniscus T,t2, the boundary condition is 



where the indices A and B refer to the fluids La and Lb, n2 is the unit 
vector of Ttt2 external to Lb, and is the mean curvature of 'Et2- 
The boundary conditions on S^i introduce an additional unknown scalar 
expressing the action of the solid wall on the fluids (sec Eq. (7), in part 1). 
The only pertinent parameter in the vicinity of the contact line is the appar- 
ent dynamic contact angle, formulated as an implicit function of the contact 
line celerity (see Section 5, in part 1); the intrinsic contact angle (see Eq. 
(10), in part 1) does no longer appear in our continuum mechanics point of 
view. Consequently, we denote simply by 6 the apparent dynamic contact 
angle and, in the following, we refer to it as the Young angle. The dynamic 
Young-Dupre equation on the contact line is 



where ly is the line friction and u the contact line celerity. The functional 
representation of v and its value were studied for a plane two-dimensional 
motion in part 1 and can be extended to an axi-symmetric motion. The 
line friction depends on the apparent dynamic contact angle (see Eq. (27) 
in part 1), and is positive. 

It is the problem of the meniscus that determines the moving boundary 
including the contact line. To this end, a complete picturing of the fluid-fluid 



pai+p^i = (l)i + Qlj, 




n2i = {Qm - QBi) + (PB - pa) n2i , 



(2) 



(72 COS 9 + ai + 1/ u = 



(3) 
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interface is necessary. In fact, the geometry of the meniscus is a consequence 
of the flow of the two fluids in the vicinity of the meniscus and this flow is a 
functional of pertinent dimensionless parameters. Without wanting to redo 
what many authors have previously done, it is useful to recall the main 
results known in the literature. 

Many papers are concerned with the shape of the meniscus near the contact 
line in capillary tubes. West, [3T], was among the first to study this problem. 
Concus, [7J, presented an analysis of the static meniscus in a right cylinder. 
Analysis using the Navier-Stokes equations and comparing results with ex- 
periments of advancing interfaces in cylinders were undertaken by Hoffman, 
[15] . Legait and Sourieau, [19], Finlow, Kota and Bose, [TT], Rame and 
Garoff, [21]. Zhou and Sheng, [32], analyzed the link between macroscopic 
behaviour of the displacement of immiscible fluids in a capillary tube and 
the microscopic parameters governing the dynamics of the moving contact 
line; Thompson and Robbins, [29], performed molecular dynamics simula- 
tions in the hydrodynamics of the contact line; Dussan, Rame and Garoff, 
[To] . Decker et al, [8], gave considerations on the contact angle measure- 
ments; Voinov, [30], proposed a thermodynamics approach to the motion of 
the contact line. In all these papers, the following dimensionless numbers 
are of significance, 

Reynolds number i?e — ^ ^ ^ ^ 



capillary number Ca = ^ ^ 
A p u^d 

Weber number We = , 

Bond number B = — — — , 
0-2 

where Ap is the density difference between the two fluids, g is the accelera- 
tion due to gravity and fiQ is the viscosity in the liquid bulk. 
A vertical capillary system treated with an asymptotic simplification was 
undertaken by Pukhnachev, [23], and Baiocci and Pukhnachev, |3j. Their 
study involves Navier-Stokes fluids paired with an imposed apparent dy- 
namic contact angle as an additional condition at the points of contact 
between the free boundary and the solid wall. The problem of the dissipa- 
tive function becoming infinite near the contact line is avoided by introduc- 
ing a slip length at the solid wall. It appears that for two-dimensional or 
axi-symmetric motions, when the capillary number tends to zero, the free 
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interface tends uniformly to the equilibrium position. 

Moreover, to minimize the effect of gravity on the meniscus shape, the cap- 
illary tube radii are generally of order less than 1 mm. A quantitative 
assessment of relative effects of gravity and capillary forces can be made 
on the basis of the Bond number. In Foister, [12], it is experimentally and 
theoretically proved that the capillary effects dominate over gravity for all 
systems when B < 1. Likewise, the relative importance of inertial and capil- 
lary effects can be characterized by the Weber number. In all experimental 
systems, when We < 10~^, inertial effects did not significantly affect the 
three phase boundary motion. 

In 1980, Lowndes, [20], performed numerical calculations of the steady 
motion of the fluid meniscus in a capillary tube. He showed that when 
Re < 10-2, the meniscus formed by an incompressible Newtonian fluid can 
be determined by using the creeping flow approximation. In this case, our 
proposed non-Newtonian model has the same streamlines as the Navier- 
Stokes model, ^3j. The method used by Lowndes considers a slip length 
less than 10 Angstroms near the contact line. This length is of the same 
order as the distance from the contact line where the fluid is no longer New- 
tonian. Comparisons between the calculated and observed contact angles 
are in accordance with the experiments of Huh and Mason, p£j; conclusions 
agree with those in part 1. In the partial wetting case, when Ca < W~^, a 

Young's angle in the interval [— , — ] verifies Eq. (3). 

6 6 

In what follows, we wish to simplify the formidable mathematical problem 
and describe the dynamics of the air-liquid interface by a suitable approxi- 
mation with as few parameters as possible. The fluids are considered to be 
Newtonian except in a close (molecular) vicinity of the contact line, |13j . 
Boundary condition on the meniscus is the classical equation (2). For slow 
movements, when Re < 10"^, Ca < W~'^,We < 10~^ and B < 1, the capil- 
lary flow is axi-symmetric, and the shape of the surface Ej2 at a given time 
t is described as a spherical cap where the angle with S^i is the apparent 
dynamic contact angle. 

At 20° Celsius, for water, in c.g.s. units, p = 1,/iQ = 0.01 and for an air- 
water interface ct2 = 72. When a = 0.1, the above conditions on Re, Ca, We 
and B are largely verified with a velocity u of some centimeters per hour. 
Similar results are obtained for glycerol {p = 1.26, /Uq = 11.8 and a2 = 63) 
with a velocity u less than a meter per hour. 

In the developments of section 4 below, these conditions are assumed to be 
fulfilled. 
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4 The motion of the contact line 



Let us present some analytical formulae related to the geometry of the con- 
tact line, see fig. 1. 

The volume of the spherical cap with height h and span 2a is 

V = — (/i^ + 3a^), where h G [—a, a]. 
6 

The volume of the liquid in the cylindrical tube is constant, so 

Vo = TTo^ye + ^ih^ + 3a2) 
o 

is constant. Let Iq be such that -ko^Iq = Vq. Then, 

ye-lo = {h^ + 3a2) and tan (---) = -- . (4) 

If A = h/a where A G [—1, 1], then 

(iTa' 

Eqs. (4), (5) establish the connections between and h and between y, 
and 9. Prom y^ = z — L and taking Eq. (4) into account, we deduce 



cos 9 = —-;-r~T2\ ■ (^) 



z = lQ + L-^{h^ + 3a2) = lo + ^ - (A2 + 3). (6) 



For a given value oft and for A belonging to [—1, 1], z is a decreasing function 
of A with values in the interval [Iq + L — 2a/3, Iq + L + 2a/'i]. 
With u = dz/dt and the above expressions, the dynamic Young-Dupre Eq. 
(3) yields 

dz 2 A 

where u is also a function of z through the Young angle 9. Taking Eq. (1) 
into account, we obtain 

— = f{z,t) + K sin — with f{z,t) = crio and K = kaio . 

dt E Q, (1 + A) 

(7) 

Remarks: The surface of the spherical cap is tt (/i^+a^). The capillary energy 
of the total system per unit length of the circumference of the capillary is 
We] so, 

27raWe = 7r(72 {h^ + a^) + 2ira J a\ dz, 
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Relation (1) implies, 



= 0-2 (— + -) + / cJio dz - I kaio sin — dz, 
2a 2 Jzi 

where zj denotes the abscissa of the point I (see fig. 1). Consequently, 

We = Wo + Wl + &^ 

with 

/l2 Z 

Wo = (T2 — + aw{z - L), Wi = k(TiQea cos — , (8) 
za ea 

and C**^ = a2a/2 — k ctiq e a cos {zj /ea) is an additional constant. 
By taking relation (6) into account, we obtain 

dh 2a2 , dWo , dWi , z 

= ~Tq~, — 2 ^'^'^ -^ = -J{z,t), -—— = -kaiQsm.— . 
oz n,^ + oz oz ea 

Eq. (7) can now be written in the form 

dz dWe 

which is an equation for a unidimensional motion of a mechanism with linear 
friction where the effects of inertia (associated with the second derivative 
d?z/dt^) are absent. In subsection 4.1, we will see that the behaviour of the 
solutions of Eq. (9), when the parameter e is vanishingly small, is completely 
different from what happens when the effects of inertia are present in a 
dynamics equation. 



4.1 Asymptotic analysis of the motion of the contact Hne 

The differential equation (9) yields the motion of the contact line. The 
behaviour of the solutions when e tends to zero has already been studied 
in the literature within the framework of a problem associated with shape 
memory alloys (Abeyaratne, Chu and James, [2j). 

Note that Eq. (7) implies df/dz {z,t) = (4c72/a) (A^ - 1)/{X^ + 1)^; for 
t given and A G [—1,1] , f{z,t) is a decreasing function of z. We consider 
the case associated with an average surface tension cJio of partial wetting 
corresponding to an average Young angle 9m such that, 

CJ2 cos 9m + o"io = 0. 
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Consequently, the inequality, | crio/c2 I < Ij is verified. Since 02 is positive, 
this relation is equivalent to the inequalities 

o"2 + o"io > and 02 — oxq > 0. (10) 

The equation of the contact line motion is only valid for partial wetting, 
nevertheless, we obtain the following limit relations associated with the limit 
values of A: 

Eq. (5) expresses the connection between A and 9; analytically, for A = — 1, 
cos 6 = 1 and the Young angle 9 = 0. Eq. (6) expresses the connection 
between A and z; for A = —1, z = Iq + L + 2a/3, for which Eq. (7) implies, 

f{z,t) = -fJio - (72. 

In the same way, for A = 1, the Young angle 9 = tt and z = Iq + L — 2a/3, 
for which f{z, t) = —aio + (T2. 

Consequently, under conditions (10), for K = kaio sufficiently small, the 
differential equation (7) fulfills exactly the conditions of the asymptotic theo- 
rem presented in the Appendix. Thus, it is possible to deduce the behaviour 
of the solutions from Eq. (7) when e tends to 0. Our aim is not to discuss 
the general solution of Eq. (7) when L{t) is an arbitrary function of t. We 
only consider two significant cases encountered in experiments (Raphael and 
de Gennes [25]), namely (i) piston at rest and (ii) piston in uniform motion 
{L{t)=vo t). 

With = z — L{t) = Iq — (a/6)A(A^ + 3), let us define the function F by 

2 A 

F{ye)= f{z,t) = -2-O-2-O-10 = -cr2COs9 - aiQ. (11) 

(1 + A j 

Then, satisfies the relation 

dy^ z 
u — +uvo = F(ye) + A sm — , 
dt ea 

where ye(0) = yo and is a function of y^. Now, straightforward adapta- 
tion of the asymptotic theorem in this simple case yields the macroscopic 
behaviour of the contact line: 

When e ^ 0, ye{t) converges uniformly to y{t) € C"^[3?"'"] and satisfies the 
differential equation 

v{y) ^ = G{y) , (12) 
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where y(0) = i/q, v is now a function of y in place of and 



[F'^{y)-K'^]2 -vvQ, if y<y-, 
G{y) = { -^VQ, if y-<y<y+, (13) 

- [F2(y) - _ j,^^^ jf y+ < y . 
y^ and are constants verifying the relations 

F{y~) = K and F{y+) = -K . 

Now, we study the two main classes of the dynamical systems (Penn and 
Miller, [22]): those in which the interface is in non-equilibrium and moves 
to an equilibrium, and the others in which the advancing interface is driven 
on the solid wall of the cylinder with a constant velocity. 

4.2 (i) Piston at rest 

In this subsection, the piston is fixed {vq = and we take z = y^). In 
Eq. (8), the potential Wq is now independent of t and a convex function 
such as that shown in fig. 2. The differential equation corresponding to the 
asymptotic behaviour of Eq. (9) when e — > is expressible in terms of the 
potential Ws as follows 

aw. 

where G{y) is defined in (13). When K is sufficiently small, the driving force 
G{y) has square root singularities (see (13) with = 0). The functions Wq 
and Ws are convex, but Ws is not the limit of W^ when e ^ (as Wq 
is). When changing the scales - i.e. when e tends to zero - Ws cannot 
be regarded as the limit of the sum of separate energies associated with the 
mean energy Wq and the energy of the perturbation We — Wq. This shows 
that by changing the scale we lose the additivity property of the energy for 
the solution of the limit differential equation: Wg is not the potential energy 
limit. However, a physical interpretation of the previous limit behaviour 
can be given. When e — > 0, the potential We admits on [y^jV'^] a large 
number of local minima whose respective distances converge to zero with e. 
For any initial position yo S [y^ ,yf] the closest local minimum is reached in 
an infinite time. When e tends to zero, any initial position yo of the interval 
is located between two local minima whose gap tends to zero with 
e. On a macroscopic scale the local minima find themselves together with 
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Figure 2: Traces of the three potentials Wg, Wq, Ws when the piston is at 
rest. In suitable units, we can take the values a = 1, K = 0.005, aio = 
—0.5, (72 = 1. The curves are centered around z = Zm associated with 
6m = ArccosO.5 (coordinate z is defined up to an additive constant). The 
potentials Wq, We and Wg have been shifted up or down for better visibility 
(since potentials are defined up to an additive constant). Note that Wq 
is the limit of Wg when e tends to 0, but this is not so for Wg which is 
fundamentally different: tangent lines at the limiting points z~ and z~^ arc 
parallel to the z-axis and the function G{z) = —dWg/dz is of square root 
type at these points. 

the initial value which becomes a position of equilibrium. Note that the 
angle ^yi(O), given by the value A, drawn from Eqs. (5), (6), corresponds to 
the y~ -position. In the same way, an angle 6r{0) corresponds to the y+- 
position. In fact K = kaio is positive and, consequently, 9a{0) and 9r{0) 
are solutions of the relations 



An angle 6 € [9 ^{0) , 9 a{0)] corresponds to an initial position yo ^ [Ue iVt]- 
This is not so for yo ^ [Ue^yt]- Indeed, in this case, the second member G{y) 



Then, 



9r{0) <9m< 9a{^) . 
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of the asymptotic equation of motion (12) of the contact hne, constitutes 
a non zero attractive force toward the points y~ and The points y~ 
and yf are now reached in finite time due to the explicit convergence of the 
solutions of Eq. (12) (see Appendix for the computation of these times). 
The position y~ (resp. yf) is the position of equilibrium associated with 
yo < y~ (resp. yo > yf) to which angle 9^(0) (resp. 0i?(O)) corresponds. 
This static case indicates that the Young angle 9 is included in the interval 
[0 niQ) , 9 a{^)\- The final value of 9, denoted by 9f, depends on the initial 
position yo of y ^^nd, consequently, on the initial value 9q of 9. We obtain 
the asymptotic behaviour 

^oe (^r(0),^a(0)) =^ 9f = 9o, 

00 < eniO) ^ 9f = 9r{0) , (15) 

eo > 9a{0) =^ 9f = 9a{0) . 

4.3 (ii) Piston in uniform motion 

When vq ^ 0, it is easy to prove that for all initial conditions yo, a constant 
solution to Eq. (12) is reached in a finite time. In the Appendix, an order 
of magnitude of this time is calculated. We obtain u = vq. 

When the piston advances, vq = u > and when the piston retreats, 
Vq = u < 0. The constant solutions of Eq. (12) are given by 

F'^{y)-k^ g\q-u'^u^ = () . 

Taking into account the definition of F in relation (11), the value of the 
advancing angle 9a{u) is 



9a{u) = Arccos ' ' V k^-io + - 'j^ j . (16), 
The same arguments yield the value of the retreating angle 9r{u) 

9r{u) = Arccos ' + V^^^ ^lo + ^ ^ , (16), 

and we deduce the inequalities 

ORiu) < 9r{0) <9m< 9a{0) < 9a{u) . 
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The Young angle ^ is a function of the contact-hne celerity u in an universal 
form. In fact, the line friction u depends on 9 but v/vq^ (where vq is a 
constant friction value) belongs to the interval [1, 1.7] when 6 € [30°, 150°] 
(see part 1). A simple approximation of (16) consists to consider z/ as a 
constant (for example an average value of the line friction like 1.35 vq). 



(a) (b) 




Tl \o X To T 



Figure 3: Values of the Young angle as functions of the celerity u of the 
contact line. We plot different formal cases associated with relations (16)i 
and (16)2 and convenient values of ctiq, 02-, k, v- In panel (a), A; = and the 
Young angle 9^ is 90°. In panels (b), (c) and (d), /c / and the Young 
angle 9,n has values 90°, 120° or 60°. In cases (b), (c) and (d), we note a 
hysteresis phenemenon due to the discontinuity of 9 for u = 0. The unit of 
the u-axis depends on the fluids. 
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5 Comparison between results, experimental data 
and other models obtained in the literature 



We propose a model of hysteresis of the contact system at the interface 
hquid-fluid-sohd with the aid of which experimental results can be inter- 
preted. Moreover, we find that the behaviour matches that obtained by 
kinetic molecular arguments in the literature. 

5.1 The dynamic line tension behaviour 

The differential equation (7) is equivalent to Eq. (3) in the form 

dz 

—r = —o'2 cos u — a\ . 
dt 

Following the expression for the solid-surface energy, the emerging equation 
admits a macroscopic behaviour expressed by the differential equation (12). 
When the piston is at rest, the macroscopic contact-line motion is governed 
by the equation (see Eq. (12)): 

r.^ = ±[(a2cos0+aio)'-fe2a?o]^ , (17) 

where the sign + or — depends on the direction of the line motion and 
k cTio = max | ai — aio | is the maximum of the fluctuations in the fluid- 
solid energy with respect to its average value ai. 

The average value of the surface tension is r = cj2COS0 -|- a\Q. It has the 
dimension of a force per unit length. If Tm = —kcTio represents the value 
of r at ^ = ^a(O) (see relation (14)), when the angle 9 is close to ^a(O), we 
obtain 

= (r'-r^)^«V27;;:(r-r^)^ ■ (18) 
When the Young angles are small, on expanding to order 2, we obtain 

02-0^(0) = /?o(r-r^) , (19) 

where /9q is a suitable constant. 

The results are extendable to the case of a piston in advancing motion. Eq. 
(17) remains unchanged, but the dynamic angle of contact is such that 

{(72 cos 6 + (Tio)^ — k^CFlQ — 1^^ = . 

Noting that = — \/ (o"2 cos 9 + (Tiq)^ — u'^, the previous results are un- 
changed but I Tu \< \ Tm I- 
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When the contact Une retreats, it is easy to present similar calculations and 
to obtain analogous results. 

These results are similar to those obtained in the literature by molecular 
statistics or single defects (Ruckenstein and Dunn, [26] )• For example, Eq. 
(18) and Eq. (19) and an analysis of their consequences are also presented 
in Raphael and de Gennes, [25], Joanny and Robbins, [18j . 




5.2 Limit velocities of the contact line 

Limit velocities for advancing and retreating contact lines have been exper- 
imentally found and described in detail (Hoffmann, [15], Dussan, [9], Blake 
and Ruschak, [4j, Chen, Rame and Garoff, [6j, Decker et al, [8]). These ve- 
locities are generally outside the domain of validity of Eq. (3). Nevertheless, 
we formally extend our calculations to the case 6 S (0,7r), and the asymp- 
totic behaviour obtained in subsection 4.1 makes it possible to calculate the 
limit velocities! Of course, this extension corresponds to the fact that the 
forms of the graphs of ^ as a function of u presented in fig. 3 are similar 
to experimental graphs proposed in the literature. Indeed, the advancing 
Young angle must be smaller than vr. The celerity of the contact line is u^^ 
and the associated line friction ly.,^. Relation (16)i yields 



In the same way the retreating Young angle must be larger than 0. The 
celerity of the contact line is now uq and the associated line friction i^q. 
Relation (16)2 yields 



Notice that the velocities uq and do not have the same absolute value. 
With this crude approximation, when the line friction is chosen with a con- 
stant value ly, knowledge of (T2, ctiq = 0^1, uq and allows us to determine 
K = k dio and u. 

In the Appendix, it is proved that the representation (1) for the sur- 
face energy is only a convenient way to consider calculations with surface 
heterogeneities. The representation is extendable to any periodic function 
with variations on short intervals with respect to macroscopic sizes. Con- 
sequently, the previous results make it possible to investigate the surface 
quality and line friction by simple measurements. 






(20) 
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5.3 Connection between the dynamic contact angle, the Hne 
celerity and the line friction 

The purpose of this paragraph is to show by comparison with simple exper- 
iments that our model leads to qualitative and perhaps some quantitative 
coincidental behaviours of the contact system with experimental evidence. 
The simplest way is to consider the surface inhomogeneity given by our 
model represented by relation (1). The values of cJio, o"2 and k, allow us to 
draw the graphs of the applications given by (16)i and (16)2, (see fig. 3). 
Only the ratio (Tio/o'2 and the value of k are important. 

We present in fig. 4 experimental layouts drawn in the literature for an 
advancing motion of the contact line (Zhou and Sheng, [32]). The similarity 
of the theoretical graph and the experimental data is striking. Nevertheless, 
we must note that the experimental data in the literature is always pre- 
sented in logarithmic scales. This comparison makes it possible to obtain 
numerical values for the line friction (as an average) and the limit veloc- 
ities uq and Utt- We use the experimental results for liquid flows in contact 
with air in capillary tubes. Results do not take any explicit account of the 
inhomogeneity of the tube walls. In the model suggested by relation (1), the 
solid surface inhomogeneity is represented by the factor k. In fact, A; <C 1, 
and in relations (20) and (21) the term k"^ ct^q is neglected in comparison to 
the terms (cj2 — cio)^ and (cj2 + cio)^- 

In order to be in agreement with the experimental conditions, we suppose the 
diameter of the tube to be about 1 mm. We deduce the numerical relations 
between the celerities Ca and u for two liquids. These expressions are given 
in c.g.s. units. For water we obtain u = 7200 Ca and for glycerol u = 4.2 Ca- 
The experimental curves yield limit velocities of the contact line. In various 
measurements they represent values for Ca and range between 10~^ and 
10~^. For an air-glycerol interface, these values correspond to the limit 
celerity with a value between 0.042 cm s ~^ and 0.42 cm s and for 
an air-water interface the limit celerity lies between 72 cm s ~^ and 720 
cm s They are widely out of the range of validity of our model, but for 
a wonder, they are in good agreement with experimental results! 
For glycerol and a limit celerity of 0.42 cm s and for a static wetting 
angle Oa{0) of about 50°, we obtain a line friction of about 240 Poise. For 
water and a limit celerity of 100 cm s the line friction is about one Poise. 
Let us note that the case of water is obviously less realistic than the case of 
glycerol. These values of the line friction are of the same order of magnitude 
as those obtained with relation (26) in part 1. 
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Figure 4: Apparent dynamic contact angle plotted as a function of Ca = 
fiQu/a2- The Ca-axis is based on a logarithmic scale. Panels (a) and (b) are 
associated with the data of J. P. Stokes et al. (ref. 13 from Zhou and Sheng, 
[32]). Panels (c) and (d) are associated with the data of G.M. Femigier and 
P. Jenifer (ref. 12 from Zhou and Sheng, [32]). Solid curves are calculated 
with (16)i. The solid surfaces are assumed to be sufficiently smooth for the 
hysteresis to be small. It is easy to adjust the values of uio, (72 and u to fit 
the experimental data given by points. 

6 Conclusion 

We proposed in this paper a dynamic model of slow movements of the contact 
system between two fluids and a solid surface. Comparison between our 
results and recent experiments or behaviour inferred from statistical physics 
shows good agreement on the qualitative level and, more unexpectedly, on 
the quantitative side. The most significant innovation in this paper is the 
introduction of the notion of line friction. This term is essential to the 
construction of the model. The line friction depends on the Young angle, 
but an order of magnitude for its average value can be obtained by using 
experimental measurements in the literature. 

The Young-Dupre relation (3) takes the inhomogeneity on the solid surface 
into account. It describes the microscopic behaviour of the Young angle. 
The inhomogeneity is distributed at distances between a few tens to some 
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hundred Angstroms. This distance is that of the operating ranges of inter- 
molecular forces which command the surface energies. At the lower part of 
this scale, energies are homogenized and these distances on a macroscopic 
scale are no longer significant. The shape of the surface is prescribed on 
larger scales. 

The results are independent of the radius of the tube. Indeed, the hysteresis 
behaviour solely depends on the physico-chemical properties of the solid 
surface. The universal form of the hysteresis loop makes it possible to discuss 
the general case independently of any particular apparatus. Relations (20) 
and (21) can be written without difficulty using the inhomogeneity of the 
solid surface in a form different from Eq. (1). 

Not all the efforts on the contact line have the same effects: the one associ- 
ated with the rapid oscillations of the surface energy of the wall, (if sin (z/ea)), 
produces a work that tends to zero when the wave length of the oscillations 
tends to zero. This effort appears in the macroscopic expression of the con- 
tact line motion represented by Eq. (12). It is noteworthy that the weak 
differences between the potentials Wq and We have huge implications. This 
may seem surprising at first sight. It is due to the fact that the contact line 
is massless and consequently its motion equation is not in the same form as 
for material systems. Some authors model contact lines as lines with matter 
(Slattery, [27j): nevertheless, the inertial force associated with the mass of 
the line is generally of an order smaller than the magnitude of the force due 
to the line friction. 

Finally we note that jumps on the inhomogeneity are considered by Jansons, 
|17j . A shift factor is proposed by Hoffman, [15j . correcting the relation 
between u and 9. These considerations appear unnecessary in our model 
where movements are slow and 9 is the apparent dynamic contact angle. 
Furthermore, the experimental literature notes the influence of evaporation 
on the relaxation time for approaching the apparent dynamic contact angle 
(Penn and Miller, [22]). Similar inferences are drawn for equipments which 
are subjected to vibrations (Marmur, [21]). The relaxation time of the 
apparent dynamic contact angle is obtained by phenomenological methods 
(Hoffman, [15], Penn and Miller, [22]). It corresponds to the calculations 
carried out in the Appendix. 

Lastly, other models for high-speed motions of the contact line using a notion 
of dynamic surface energy are also considered in the literature: Shikhmurzaev, 
[28] , Blake, Bracke and Shikhmurzaev, [5]. Their study is outside the scope 
of our slow movement model. 
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7 Appendix: Proof of the fundamental theorem 

Consider the differential equation 



^(^) 37 = /(^' t) + K sin - 
at e 



(Al) 



(we present the same differential equation as Eq. (7), but / may be a more 
general function than listed in (7), and o = 1 corresponds to the choice of a 
convenient unit of length) . 

The following hypotheses are assumed: 

• K and e (e <C 1) are two strictly positive constants; 

• / is continuously differentiable for any z G 3?, and for any t G 

df/dz < ; 

• There exist z~ and z"*" belonging to C^(3f?) such that for any t S 



f{z-,t) = K and f{z+,t) = -K- 
is a strictly positive continuous function of z. 



iA2) 



Remarks: The fact that for t fixed, f{z,t) is a decreasing function on 
[0, +00), implies z~ < z'^. Differential equation (Al) yields a single solution 
Zs{t) defined in 3ft+ with Z£{0) = zq, (the differentiable equation fulfills the 
conditions of uniqueness of the Cauchy problem (Hartman, |14j)). 
We obtain the fundamental result that gives the behaviour of ^^(i) when e 
tends to from the following theorem 

Asymptotic theorem 

When e tends to 0, Zc{t) converges uniformly to z{t) belonging to C^[3ft+] 
and satisfying the differential equation 



where z{0) = zq and 



dz 


= 9 


z,t) , 




if 


z < z'{t) , 





if 


z-{t) < z< z+(t) 




if 


z+{t) < z. 



(AS) 
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The theorem, proved for ^{z) = constant, is extended without difficulty for 
^{z) befonging to [a, /3] where 0<Q!</3<cxd[1. 

Let us give an elementary proof of the fundamental theorem. For a 
complete demonstration using Young measures, we refer to Abeyaratne, Chu 
and James, [1], [2]. We just consider the case i' = constant (z^ = 1 with 
a convenient unit) and / independent of t. For the variable z belonging 
to a compact interval, the segment [—1,1] is considered. The hypothesis 
df /dz < corresponds to strict convexity of Wq (j!)Wo/dz = — /(z)^ . For a 

given e, consider zq £ [—1, 1] as the initial value of the solution Zs{t) of Eq. 
(Al). The critical points of Eq. (Al) are the roots of 

f{z)+Ksm^ = 0. 

The roots belong to the interval [z~ , z~^]. Let us consider the case for which 
Zq is smaller than z~; when T is sufficiently small, it is the same for Ze(t) 
with t S [0,T]. Let us calculate the limit values t according to position z 
when e tends to zero. For e ^ 1, the function sin(z/e) oscillates rapidly 
between -1 and 1. Let zq.e be the value immediately above or equal to zq 
such that sin(zo,£/e) = — 1 and let zi^e be the value immediately below or 
equal to Zs{t) such that sm{zi^i;/e) = 1. Then, 

lim / ^ = t . 



2; 

20, e f{z) + Ksin- 



Let us divide the interval [zq^s, zi^^] in intervals of length evr. Then, 



/ 



dz _ /■2o,.+(p+l)e^ 

2o,e f{z)+Ksui- J ZQ^s+peiT /(z) + ivTsin- 

where = (zi,^ - zo,e)/(£7r) 

At each interval [zq^^ + pevr, zo,^ + {p+ l)£'?r] the change of variables s 
sin - yields 

£ 

dz f^^ eds 



£ 



^If £, ~ P i'{(;)d<^, ^ is a strictly increasing function of z corresponding to a bounded 
change of length scale depending on the considered point and fi{(,,t) = f{z,t) is a de- 
creasing function of 5. 
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vr vr 

where A = earcsins and A G [epir — e—,ep7T + e—]. Because / is continuous, 
this expression has the same hmit when e tends to zero as 

ds 



p=0 

or 



/ £7r , , 

T^t'n J-i {f{ep7r-e^) + Ks)VT^ 



1^' ds 



en 



J-1 {f{epn + e^) + Ks)Vl^' 
which are two sums of Darboux integrals 

ds 



i ( f 

Jzo,e ^ \J-1 



{f{z)+Ks)^/l^ 



dz . 



When £ tends to zero, this expression converges to 
1 / ds 



dz = t . 



{f{z)+Ks)Vl^, 
Because f{z)'^ > K"^, the change of variables s = sina yields 



Li {f{z) + Ks)^/l^ ^f{zf-K^ 

and finally, 



dz 



which implies 

In the same way, for zq above z'^ we obtain, 

i = -V/(^)2-K2. 

When Zq belongs to the interval [z~ , z^], it is easily seen that there exists a 
critical point 1^ of Eq. (Al) such that \z^ — zq \< be where 6 is a positive 
constant depending only on K and /. The solution of Eq. (Al) tends in 
a monotonous way towards 1^. Consequently, for t belonging to [0,+cx)[, 
Z£{t) converges uniformly to zq when e tends to zero. The macroscopic law 
associated to zq belonging to [z~ , z~^] is i = 0. 
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General oscillations 

We consider the superimposed effect of an arbitrary smooth periodic func- 
tion P[z/e) in place of sin(z/e). We can assume without loss of generality 
that P{z/e) has a zero average, (in other cases, we add a constant to /). 
The amplitude of P gives the placement of the flat region as in figure 2: 



With the same hypotheses as in the theorem we define z (t) and z^{t) such 
that 

f{z'{t),t) = -minP and f{z^{t),t) = -maxP ; 
so the flat region remains [z~ [t), z^ {t)]. 

In the special case when / is independent of t, the calculation is developed 
in the same way as previously. The results are unaltered with — minP in 
place of K and — maxP in place of —K. 

The complete proof of this extension is given in Abeyaratne, Chu and James, 



Relaxation time associated with Eq. (A3) 

We consider the case when g is explicitly independent of t and v is constant. 
For zq < z~ , let 



For zq near z , f{u) + K ^ 2 K and f{u) — K ^ {u — z )f'{z ). Then, 
/(u)2 - ^ 2K{u - z-)f'{z-) and 



The x-value yields the magnitude of the relaxation time necessary to obtain 
the final position of the contact line. 
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